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Abstract 

A simple lattice model inducing a gauge theory is considered. The model de- 
scribes an interaction of a gauge field to an N x N complex matrix scalar field trans- 
forming as a field in the fundamental representation. In contrast to the Kazakov- 
Migdal model the model contains only the linear interaction between scalar and 
gauge lattice fields. This model does not suffer from extra local C7(l) symmetries. 

In an approximation of a translation invariant master field the large N limit of 
the model is investigated. At large N the gauge fields can be integrated out yielding 
an effective theory describing an interaction of eigenvalues of the master field. The 
reduced model exhibits phase transitions at the points /3 c - r and (5 cr and the region 
(Per, Par) separates the strong and weak regions of the model. 

To study the behaviour of the model at large N in more systematic way the 
quenched momentum prescription with constraints for treating the large N limit of 
gauge theories is used. With the help of the technique of orthogonal polynomials 
nonlinear equations describing the large N limit of the reduced model with quenching 
are presented. 
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1 Introduction 



This paper is inspired by recent studies of the Kazakov Migdal (KM) model [T| - [14|, which 
was suggested as a possible model of hadrons. The KM model contains a scalar field $ 
in the adjoint representation of the SU (N) group interacting with the usual lattice gauge 
field. The main point which allows to study the large N limit of the model analytically, 
is the absence of the kinetic term for the gauge field. 

There are several problems with the KM model. A relation of this model with QCD 
and, in particular, the property of asymptotic freedom still should be clarified. The KM 



model has an extra local gauge symmetry 0, Oj, which should be broken if we want to 



recover QCD from this lattice model. There is also a more technical problem, concerning 
the method of investigation of the large N limit in the model. Most of previous consid- 
erations used the approximation of the constant master field. To clarify this assumption 
the large N limit of the KM model has been investigated by an other method, namely, 



by the quenched momentum prescription [15| - [Eq] for study the large N limit of gauge 



theories |24|. It has been shown that using the quenched procedure with constraints one 
obtains the one side model with the quartic interaction, that makes the model not exactly 
solvable. 

In this paper we are going to construct a modified model which does not suffer from 
extra gauge symmetry and admits more simple analytical treating at large N. 

Several modifications of the KM model have been proposed recently. Migdal [ |l2|j , 
and Khokhlachev and Makeenko proposed to induce QCD by fermions in the adjoint 
representation of the gauge group SU(N) with various types of lattice fermions. They 
discussed also adjoint scalar and fermionic models at large number of flavours Nf. 

We will consider a simple lattice gauge model in which an interaction of gauge fields 
is induced by a scalar field, being a.n N x N complex matrix transforming as a field in the 
fundamental representation. In contrast to the KM model the modified model contains 
only linear interaction between scalar and gauge lattice fields. This model does not suffer 
from extra local U(l) symmetries. To estimate the behaviour of this model we perform 
an integration over lattice gauge fields and then use a translation invariant master field 
to calculate the remaining integral over scalar field. The integration over gauge fields at 
large N is performed by using the Brezin-Gross results [26 of study one link gauge model 
in an external matrix field. The phase transition which takes place at large N for one link 
gauge model in an external matrix field [p7|, |26[| gives rise to a non-trivial phase structure 
of our model. The reduced model exhibits phase transitions at the points /3 c - r and Pgr and 
the region (/? c - r , /?„■) separates the strong and weak coupling regimes of the model. 

To study the behaviour of the model at large N in a more systematic way the quenched 
momentum prescription with constraints for treating the large N limit of gauge theories is 
used. This prescription leads to a quadratic dependence of the reduced action on unitary 
matrix V instead of a quartic dependence for the KM model. This quadratic dependence 
permits to use the technique of orthogonal polynomials [p9| [l[ and to obtain nonlinear 
equations describing the large N limit of the reduced model with quenching. These equa- 
tions are similar to equations describing the KM reduced model without quenching. 
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2 The Action of the Model 



Let us consider a simple lattice model which describes an interaction between N x N 
matrix scalar field $ and gauge field and U^. The partition function is given by 

Z= [Y[dU„(x)Yld$(x)exp{S(<Z>(x),U^x))}. (1) 

X,[l X 

Where $ is an arbitrary N x N matrix and the scalar field <&(x) defined on the sides 
of a d-dimensional rectangular lattice and is an unitary matrix defined on the links; 
<i$ = Ili.j d(pij an d dUfj, is the Haar measure on the group of unitary matrices. The action 
has the following form 

S($,U) = --m 2 N tr $$t +(3NY, tr [$(x)[/ M (x)$ t (x + /i) + + ^)[/J(x)$t( x )]. (2) 

As in the case of the KM model in (Q) there is no a kinetic term for the gauge field C/ M . 
The action (0) is invariant under gauge transformations 

U^x) -> fi(x)^(x)fi f (x + /i), $(x) -»• ^(^(z), $(x) f -> n(x)$ f (x) (3) 

Note that one can consider the matrix field 0^- as a field with upper and down indices 
4>ij = (p{ and regard the upper index as the usual colour index and the down as the flavour 
index. So, with the respect to the down indeces the scalar field transforms as the usual 
matter field in the fundamental representation, and the action is the sum of N copies of 
matter actions in the fundamental representation. In accordance with this analogy the 
action is invariant under the global U(N) symmetry 

$(x) -> G*$(x), -> &(x)G, (4) 

G belongs to U (N) . Note also that in contrast to the KM model the model (|2|) does not 
have extra local U(l) symmetries. 

Calculating gaussian integral over fields $ and one rewrites the model @ as a pure 
gauge theory with the following reduced action 

S md = NU^)-^hrU(T) (5) 

Notice that the expansion ([5]) starts with the standard Wilson action. Recall that the 
analogous expansion for the KM model starts with the square of the Wilson action. To get 
the Wilson action instead of the action (|5|) one has to find a mechanism of suppression of 
long loops in the effective action (H). For example, one can use the following trick, which 



was also used by Khokhlachev and Makeenko JT3[ to the KM model. Let attache an extra 
index to the field $, i.e. introduce the field , / = 1, ...rif and consider mo — > 00 and 
n/ — ► 00 so that rriQ ~ (n/) 1 / 4 . Then only the single-plaquette survives in (|5|), the longer 
loops are suppressed at least as nj 1 ^ 2 . 

A naive continuous version of our model corresponds to iV copies of the gauged scalar 
fields in the fundamental representation. These N copies can be interpret as fields with 
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different ffavour and N c = Nf. The corresponding Lagrangian has no an explicit Yang- 
Mills term 

C = — (<9$ - i$A u )(d$*) + iAJ& + V($&). (6) 

Here there is not the problem which arises in the naive continuous version of the KM 
model, where a commutator [A^, $] presents and therefore we miss the diagonal degrees 
of freedom of the gauge field. 



3 Reduced Model 

In this section we estimate the partition function flU) in the large N limit using the 



quenched momentum prescription WB - p3[. This prescription for treating the limit 



of infinite N of matrix field theories has been obtained more than ten years ego by Eguchi 
and Kawai [T5|], Bhanot, Heller and Neuberger [ITH], Parisi fT7|, and Gross and Kitazawa 



18fl (see also |L9| - [p3| ). Generally speaking, according to this prescription one gets 



the vacuum energy in the large N limit by integration the free energy of the quenched 
reduced model over the quenching parameter p. The reduced model is described by a 
reduced action and to get the reduced action one should replace a matrix field $(x) with 
D(x)QD^(x), where D(x) = e tp ^ XtJ -, and p^ is the diagonal matrix with matrix elements 
p^, i = 1, ...N. (quenching parameters). It has been shown that quenched theory produces 
the standard Feynman diagrams for invariant Green functions in all orders in perturbation 
theory. 

Performing the reduction procedure $(x + /i) = D^(x)Dj 1 and a change of variables 
U — > U' = U^D^ we left with an action 

S($,U,D) = --m 2 Ntr<£>& + tr [QU^D^ + QU^D^ (7) 

2 ti>0 

Integrating over $ and $t one gets the sum over T of terms m~ 21 ^ trU (T) tr D (T)). 
Since all D are diagonal matrices for the closed contour V we have D(T) = I, and, 
therefore, we can neglect the last multipliers in both terms in equation (|7|) and we left 
with the following reduced action 

S red ($,U) = ~m 2 NtT$& + PNY, tr^U^x)^ + ^U^) (8) 

2 fM>0 

An other way of obtaining the reduced action (^j) consists in application of the reduced 
procedure directly in the induced action ([5]). Performing the reduction procedure for the 
gauge field U^x+u) = D u U tl (x)Dl we get the sum of terms U(T), this answer follows also 
from (|8]). The same action gives the approximation of the translation invariant master 
field. 

Hence, we have to estimate at large N the following free energy 

/ffb — 
d<$>dU lM exp((-—Ntv$& + /3NY^ tr + C/t)$ f )) (9) 

2 /i>0 

where f3 = a d /3, rh 2 = a d m 2 . 
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Note that in (|9]) $ is an arbitrary complex NxN matrix. Any matrix $ may be written 
in the form $ = PV with P hermitian and V unitary. P may be diagonalized by a unitary 
transformation, P = SAS^. If we perform the change of variables U — > S^VUV^S in the 
integral 

J($,$t)= [dU lt exp{PN y EK*(U lt + Ut)$) (10) 

it is then clear that the integral (|T0|) depends only upon A. Up to a normalization factor 
the measure d$ has the following form d 2N $ = A 2 (A) UidXidSdV. The integral (JT^) 
defines the well known object, this is the partition function of the Brezin-Gross model 
which describes one link gauge field in the external matrix source [^, ^6|. The form of 
this partition function in large N limit depends on the magnitude of external source. In 
our case this parameter is 

S=itr($$t)-i = i.yu|-a (11) 
Nf3 y 1 Nf3 i ' v ; 

where <pi = (A)&. The weak coupling regime corresponds to S < 2, and the strong 
coupling regime to S > 2. The answer for both regimes can be written in the following 
form p6 



4 = |E(^l0| 4 + C ) 1/2 -^Elogp 2 |0,| 4 + c) 1/2 + (/5 2 |0/ + c) 1/2 ]-(c + 3/4) (12) 

The constant c is different for different regimes. For the weak coupling regime c = p5 
and for the strong coupling regime the constant is determined from the equation 

i£(/? 2 |0*| 4 + c)- 1/2 = 2 (13) 
The next step is the calculation the integral in respect of 0's 

F = In J n(& - <f>j) 2 II d^dc[6(S - 2)8(c) + 6(2 - S)5(c - ~ £ M"" 2 )]' ( 14 ) 



i,3 



exp[-^iVE # + ^d^WW + c) 1/2 " 
d - £ logf^l^l 4 + c) 1/2 + (^ 2 |0,| 4 + of] - (c + 3/4)] 

i,3 

3.1 Week Coupling Reduced Model 

One can roughly estimate the behaviour of this model in week coupling regime ignoring 
the contribution from the second term in the square brackets in the first line of ( |H| ) and 
substituting in (|14]) an average of 

<0 2 >=ltr$ 2 = i. (15) 
J JS fj, 

instead of the sum over j. In this case we get the usual 1-matrix problem in the logarithmic 
potential 

V = — ^Ntr$ 2 + -tr log (1 + /i$ 2 ) (16) 
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where = 4d(3—m 2 . The parameter li should be defined from the selfconsistent equation 
(|T5|). Taking into account only quadratic terms in ( [TED we get the following selfconsistent 
equation 

d\x — ml = //, (17) 

which has a solution for d > 1 only if m\ > 0, i.e f3 cr = rh 2 /Ad is a critical point (the 
symbol cr means that this critical point is obtained starting from week coupling regime). 
There is also a restriction on d < 5/4 coming from the requirement of applicability of the 
used approximation. According to this requirement the region |0| < must be inside 

the region where the potential (|T6|) increases, i.e. 4//i < djm\ — 

However there is a difficulty which makes the application of the quadratic potential 
V($) in the weak coupling regime non-consistent. It is caused by a breaking of the 
equality S < 2. This gives a raison to add a quartic term to the potential 

V($) = — -iVtr$ 2 + #iVtr$ 4 (18) 

Here m 2 < Adf3, and, therefore, m 2 , is positive, i.e. we have the Goldstone potential; 
g = a d g. In this case equation (|23| ) gives li = In the adopted approximation 

the condition S < 2 means > 2, that can be achieved for j3 being large enough; 

9eff _ 

g e jf — g — This approximation has also the following restriction > 8 ^£^ ■ Note 
that all the numerical factors are obtained in the approximation of the quadratic potential 
in the vicinity of the minimum of the Goldstone potential. 

One should note that the application of the reduction procedure itself for the weak cou- 
pling regime is not quit correct since in this region as it was mentioned in the Introduction 



one has to use the reduction with constrains (see [23]). However the simple calculations 
presented above show that we have to deal with the effective Goldstone potential and this 
gives an indication that the structure of the model with quenching momentum may be 



similar to the phase structure of spin glasses [19 



3.2 Strong Coupling Reduced Model 



Now we are going to estimate the free energy (plD in the strong coupling regime. In this 
case the second term in the first line of ([24]) gives a main contribution, i.e. we have to 
deal with (|22"), where c satisfies (p3|). Taking the following approximation for (p2j) 



F=-Ng eff [3 2 Y, ( t> 4 + const > 9e ff = 2 ^f^ (19) 
we get the one matrix problem in quartic potential with a wrong sign. A phase transition 

1 /5 m 2 



in this model occurs at /5 c - r = lyf^) the point where the function of distribution of 
eigenvalues lose the positivity (here cr means that this critical point is obtained starting 
from strong coupling regime). It is interesting to compare the strong coupling regime of 
(||) with the strong coupling regime of the KM model. We will do it in the forthcoming 
publication. 

The above estimations advocate for the following picture. There are two phase in 
the theory (|), weak and strong phases. They are separated by the intermediate region 
(Per, Per) ■ The presented estimations are rather rough and in principle they can be made 
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more precise if one does not restrict by the quartic potential and deal with logarithmic 
potential. However we believe that the qualitative picture of the phase structure of the 
model (|) is correctly reproduced. The condition /3 c - r < j3cr_ gives a restriction on d. 



4 Reduced Model with Quenching 

In the case of gauge theories the quenching of the momentum must be accompanied by 
a constraint on the eigenvalues of the covariant derivative. Without this constraint one 
gets naive reduced model without quenching. In the case of the Wilson gauge theory the 
naive reduce model describes correctly the theory only in the strong coupling regime. The 
suitable constraint restricts the eigenvalues of U^D^ to be equal to _D M . This means that 
we have to take the measure dfi(U) in @ to be 

d f i(U) = I[dU li C(U li ,D li ) (20) 

where dU^ is the Haar measure on SU(N) and 

C(U„ D,) = Y[[ dV, 6(U, - VrDrfDl) (21) 
i" J 

If we integrate out the {7 M , setting U ^ = V^D^V^D^ we obtain 

exp(F(p)) = (22) 

( d$l[dV^exp{a d (~m 2 Ntr$& + tr [QV^D^&Dl + W^lV^ D „)]} 

To find the vacuum energy at the large N limit one should integrate F(p) over dp p. In 
fact we will use the formula 

F= / dv(p)F(p), (23) 



where 

dM = U-^ d eM~pW), (24) 

So, the reduced action with quenching for the model (0) contains the quadratic depen- 
dence from the unitary matrix U. Note that the reduced action with quenching for the 
KM model has the quartic dependence from the unitary matrix, that makes the model 
not exactly solvable. 

The first step in the calculation of (|22D is the evaluation of the integral over 's. It is 
interesting to note that just for this case the technics developed by Kazakov and Migdal 
fl], 0, and Gross JTT|] permits to study analytically the problem in question. 



Using the two-matrix representation of the Itzykson-Zuber (IZ) integral and the 
orthogonal polynomial method |29 the following integral over 's 

exp/({0 2 },{pi}) = J dV exp(iV tr $ (V(D + D*)V*&)) (25) 

at large N can be represented as 

exp/({0 2 },fe}) =N 2 /■ 1 dt(l-t)ln/(t), (26) 

Jo 
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where 



/(iH _ t + / «M> (27) 



y J 2m 1 — zy/(t) 
tfe K'( p(z),t) 



. 1 / dz \ 



and 



v fc) = xf EM*" 4>l),V 2 {y) = 2 cos a H ), (28) 

b b 

We left only linear dependence on D and in the action in (p5[) , since the dropped 
D's do not contribute to closed loops. Note that in spite of this reason such step needs 
more arguments. 

The next step is the evaluation of the integral over $ and p^. Using the replica trick 



17, 19 



ln/#A 2 (0)exp{/({0 2 },fe})} = 

exp{/({(<n 2 },fe})}-ll (29) 




one can perform first the integration over p^s and then use the saddle-point approximation 
for integral over <ft 

2 V = 2N-^<f)? - d^— J(U a }) (30) 

where 

l ni 7(0 Q ) = J ^expl({0n 2 },fe}), (31) 



a is the replica index. Equation (|30|) is more complicated then the corresponding equa- 
tion for the reduced KM model and at present we do not know how to solve it. However, 
an obvious advantage of the model (|l|) is that the IZ integral completely describes the 
behaviour of the model, while this integral describes the KM model only in the approxi- 
mation without quenching, and even exhaustive knowledge of the result of integration of 
the IZ integral over eigenvalues of $ does not provide information about the behaviour of 
the KM model in all regions of coupling constant. 



5 Concluding Remarks 

A lattice gauge model with simple linear interaction between gauge and scalar fields was 
considered. This model does not suffer from local U(l) symmetries and it contains the 
Wilson action as the first term in the expansion of the effective action. So, it seems the 
model can be considered as a rather suitable approximation to QCD. It was shown that 
the model in the large N limit admits the analytical investigation and has the rich phase 
structure. There are many open problems with the model. The most important question 
here, as well as for the KM model, is a clarification of the status of asymptotic freedom. 

There are also some more technical problems which must be trieted out. First, it 
would be interesting to study the reduced model in an external random field. The action 
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of this model is given by the expansion (|7]). Note that the role of the random field D 
consists in the suppression of contributions coming from open loops. A similar suppression 
may be ensured by abelian random fields exp(i</? M ). Such model will be a subject of our 
forthcoming paper RT[. Another technical questions, as it has been mentioned in Section3, 



are related with more thorough estimation of the free action (PI). 

In concluding let us note once again that though there are many open questions with 
induced QCD, this subject seems to be very promising and undoubtedly is worthy of 
further developments. 
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